Anomalous Aharonov-Bohm conductance oscillations from topological insulator 

surface states 
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We study transport properties of a topological insulator nanowire when a magnetic field is applied 
along its length. We predict that with strong surface disorder, a characteristic signature of the band 
topology is revealed in Aharonov Bohm (AB) oscillations of the conductance. These oscillations have 
a component with anomalous period $o = hc/e, and with conductance maxima at odd multiples of 
|$o, i.e. when the AB phase for surface electrons is n. This is intimately connected to the band 
topology and a surface curvature induced Berry phase, special to topological insulator surfaces. We 
discuss similarities and differences from recent experiments on Bi2Se3 nanoribbons, and optimal 
conditions for observing this effect. 

PACS numbers: 



There has been much recent interest in Topological 
Insulators (TIs) , three dimensional solids which are insu- 
lating in the bulk but display protected metallic surface 
states (see £Q for reviews) , in the presence of time rever- 
sal (T) invariance. Surface sensitive experiments such 
as ARPES and STM pQ , have confirmed the existence of 
this exotic surface metal. However, so far there have been 
no transport experiments verifying the topological nature 
of the surface states. Besides the fact that bulk impu- 
rity bands may contribute significantly to conductivity, 
the key transport property is the absence of localization 
when T invariance is present. This has been hard to con- 
vert into a clear-cut experimental test. Here we discuss 
a topological feature of TI surface states that can lead to 
a transport signature. 

Consider a wire of topological insulator, with magnetic 
flux applied along its length. The surface states can be 
considered a collection of one dimensional modes, that 
come in pairs moving up and down the wire. Time rever- 
sal symmetry is present at zero flux, and also, approxi- 
mately for the surface states, when the surface encloses 
an integer multiple of |$o (= hc/2e) flux quanta. How- 
ever, there is an important difference between even and 
odd multiples of |<&o-flux. Odd multiples of |$o-Aux 
leads to a 7r Aharonov-Bohm phase for surface electrons, 
resulting in an odd number of pairs of modes. This one 
dimensional state is topologically protected and cannot 
be localized with T symmetry[2]. In contrast, even multi- 
ples of 5 •I'o-flux leads to an even number of modes which 
are not protected. Thus, with sufficiently strong disorder, 
even flux leads to a fully localized state for a long wire 
while odd flux leads to a metallic state whose conduc- 
tance approaches 2e 2 /h, under ideal conditions. Inter- 
estingly, as discussed below, a crucial component of this 
even-odd effect is a surface curvature induced 7r Berry's 
phase. Throughout we assume low temperatures so the 
thermal dephasing length exceeds the sample dimensions. 

Note, this oscillatory dependence has a flux hc/e pe- 



riod, in contrast to Aharonov- Altshuler-Spivak (AAS) os- 
cillations which have period hc/2e. There has been much 
discussion on the question of hc/e vs hc/2e oscillations, 
in mesoscopic rings [3J ^ and cylinders E]. For rings, 
both periods are observed, the first period is attributed 
to AB interference of single electrons, and the second 
to AAS oscillation arising from weak localization effects. 
However in metallic cylinders, only the hc/2e period has 
been experimentally reported [5]. The hc/e period has 
been theoretically predicted, [H [6], but always occurs 
with random sign, i.e. they can peak at either even or 
odd multiples of |3>o- Therefore an ensemble average 
tends to wash out this effect, which, according to Ref. 
OH], is why the hc/2e effect is more commonly observed. 

In contrast, the hc/e oscillation described here are 
unique in having maxima always at odd-integer multiples 
of ^$0! and only occur in strong topological insulators. 

A recent experiment on topological insulator Bi2Se3 
nanowires has indeed reported such an anomalous hc/e 
flux period 7 J . However, there is a crucial difference from 
the effect described above - the conductivity is found to 
be minimum at the locations of the predicted maxima. 
Hence Ref. [7] is presumably observing different physics, 
explaining which is an interesting open question. The 
regime described in this paper is best accessed by going 
to strong disorder on the surface, or by enhancing the 
one dimensional nature of the system, eg. by considering 
narrower wires. We believe this regime could well be 
accessed by future experiments of a similar nature. 

We first describe the physical ingredients that give rise 
to this anomalous AB effect, such as curvature induced 
Berry phase, in clean systems. Subsequently we report 
the result of numerical experiments on disordered cylin- 
ders of topological insulators, realized in a three dimen- 
sional lattice model. Strong disorder is confirmed to ex- 
pose this anomalous AB effect. 

A single Dirac cone is the simplest model of TI sur- 
face states, which naturally invites comparison with 



graphene, with a pair of Dirac nodes centered at mo- 
menta ±K , for each spin projection. If scattering be- 
tween the graphene's nodes are neglected, one might 
conclude that the topological insulator surface is simply 
'l/4th' of graphene. Here we point out a topological ef- 
fect that is special to TI surfaces, connected to the fact 
that the surface Dirac fermions are sensitive to spatial 
topology, in the same way as Dirac particles in a curved 
two dimensional space. There is no analog of this for 
graphene, even when it is rolled up into curved structures 
such as nanotubes and buckyballs. The root of this dif- 
ference can be traced back to the physical interpretation 
of the Dirac matrices. While in graphene, Dirac matrices 
act on an internal psuedospin space, in topological insu- 
lators they act on physical spin, which are locked to the 
surface orientation. Hence surface curvature introduces 
new effects in TIs as described in Ref I8J. In particular, 
consider an electron circling the cylindrical surface of a 
topological insulator. Due to the locking of spin to the 
surface orientation, an additional Berry's phase of n is 
acquired during such a revolution. The surface modes of 
a cylinder appear in pairs moving up and down the axis. 
A consequence of the Berry phase is that there are an 
even number of these pairs. If modes are labeled by an- 
gular momenta, kg, these are quantized to half integers 
because of the Berry phase. Thus there is no unpaired 
low energy Dirac mode at kg = 0. In contrast, in carbon 
nanotubes, the kg = mode is present and responsible 
for the metallicity of eg. armchair nanotubes. Now, if an 
additional Aharonov-Bohm flux of ir threads the cylinder 
of topological insulator, canceling the curvature Berry's 
phase, one reverts to the regular quantization i.e. kg — 
is allowed, and hence an odd number of one dimensional 
mode pairs are present (see Fig. [Ta,b). 
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FIG. 1: (a) proposed geometry, with (weak) magnetic field 
applied along the wire axis, with flux p = 27r"I>/$o enclosed. 
(b) Quantization of transverse momenta kg for a cylinder, as 
a function of applied flux. The shaded area boundary in the 
Fermi surface. For <p = 0, an even number of ID mode pairs 
occurs; but with ip — n flux, an odd number of pairs, pro- 
tected by time reversal symmetry, is expected. Demonstra- 
tion in a microscopic model: spectrum of a cuboid of clean 
topological insulator, with cross section 10x10 unit cells (ge- 
ometry shown in inset), as a function of momentum k z along 
its length, (c) Spectrum for p = 0, each band shown is dou- 
bly degenerate (d) For flux p> = ir, all modes except the linear 
one are doubly degenerate; hence an odd number of ID mode 
pairs is present. 



Surface Dirac Theory: Consider generalizing the Dirac 
Hamiltonian for a flat surface perpendicular to the z di- 
rection, H = —ivph(p x a x + pyCTy) to the case when the 
surface is curved. We utilize the fact that the surface is 
embedded in three dimensional space, so t(x , x 2 ) is the 
three vector defining the surface location, as the coordi- 
nates x l [i — 1, 2) are varied. Then, e^ = dr/dx l are 
tangent vectors. Define conjugate tangent vectors e J , via 
e* • e^ — Sij, the Kronecker Delta function. Naively, one 
might guess that the Dirac equation on this curved sur- 
face is just: Hi — ~iv Fh~{a l d\ + o?d-i), where di = d x i 
and a 1 — e l ■ a are the Pauli matrices along the tan- 
gent vectors. However, the actual form is a little more 
involved, and is most compactly written if we assume x l 
are normal (also called geodesic) coordinates. In these 
coordinates, (i) the e* are orthonormal e l • e J = <5y; and 
(ii) the derivatives: de l /dx l , are along the surface nor- 
mal e 3 (= e 1 x e 2 /^ 1 x e 2 |). The coefficients of propor- 
tionality are the principal curvatures de l /dx J = — -^e 3 . 
Such a set of coordinates can always be found locally. In 
these coordinates it is readily established that there is an 



additional term: 
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H D = -iv F h{a 1 d 1 +a 2 d 2 ) + i^f- (— + — )/3 (1) 

where (3 = e 3 • a . A shortcut to obtaining it is requiring 
that the Hamiltonian H]j be Hermitian and also anticom- 
mute with the local /? matrix (note, the matrices a 1 , ft 
are space dependent). 

A more formal derivation of the same result proceeds 
from the general form for a Dirac equation in curved 
space: Ho = —ihvFi^Di + a 2 !^), where Di is the co- 
variant derivative along a pair of general coordinates x l , 
defined as Di = di + Ti, where Ti is the spin connection, 
which here is given by the equation: 



[T i ,a j ] = d i a j +T: 



tk L 



(2) 



where a k 

is assumed. The T J ik are the Christoffel symbols. This 

can be solved explicitly to give the following simple ex- 



e • <t, and summation over repeated indices 
The Tl k 



pression: 



r 4 = --fib? 



(3) 



The Christoffel symbols arc defined from the derivatives 
of the tangent vectors, projected into the tangent plane: 
[c^e- 7 ],, = —T J ik e k (indices take on values £ {1, 2})[5]. 
Now, we switch to normal coordinates. Clearly, the 
Christoffel symbols vanish since the derivatives of the 
tangent vector are now normal to the surface. Thus from 
Eqn 2 the spin connection is given by I\2 = ± 



2_Ri 



whicn leads back to Eqn. [T] An alternate elegant formal- 
ism is developed in Ref. [TU], for quantized Hall states. 

Now, let us specialize to a cylindrical surface such that 
x 1 = z along the cylinder axis and x 2 — R9, where R is 
the radius and 9 is the angle around the cylinder. Now, 
Ri = go and R 2 = R, in Eqn. [11 The unitary transfor- 
mation U = e lcrz8 / 2 , transforms that into the canonical 
form H' D = —ihvp(<7 z d z + Uydg/R). However, since the 
unitary transformation changes sign 9 — > 9 + 2w, the 
wavefunctions for the new Hamiltonian satisfy antiperi- 
odic boundary conditions on circling the cylinder. There- 
fore, only angular momenta h(m+l/2) are allowed, where 
m is integer. Hence, the zero angular momentum is ab- 
sent, and there are an even number of one dimensional 
modes pairs. Now threading an additional it flux, the pe- 
riodic boundary conditions are restored, and the parity 
of the mode pairs is reversed; see [9] for a more general 
argument. Although the cylinder has vanishing Gaussian 
curvature, a nonzero spin connection leads to the Berry's 
phase of n. This topological property is also ultimately 
responsible for metallic dislocation lines [U [11] . 

Microscopic Model: We now demonstrate this effect for 
a lattice model of a strong topological insulator (which 
is more general than the Dirac approximation). We use 
the model of Fu-Kane-Mele [TJ] on the diamond lattice 
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(Alk x d, 2 fc ) • Wc fcff '(4) 



Parameters are chosen to give a strong topological in- 
sulator [5] with bulk gap A = 2t. A long cuboid with 
cross section L x L is taken along the weak index di- 
rection of this model, and surface states are labeled by 
momenta k z along the long axis. A uniform magnetic 
flux $ is introduced uniformly through the cross section, 
denoted in units of the flux quantum: ip = 27r$/$o- The 
surface spectrum is shown in Fig. [T] All modes are dou- 
bly degenerate except the linearly dispersing mode in[T]d. 
Thus, even for small sizes L — 10, the even-odd mode ef- 
fect and the gap closing at flux ip = n is apparent. While 
the breaking of time reversal symmetry at this flux im- 
plies there is always a gap, this is seen to be very small, 
and the curves appear 2ir periodic, so time reversal is ap- 
proximately a good symmetry at these flux values. Thus, 
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FIG. 2: (a), (b) and (c): Localization length A in the quasi 
one dimensional geometry, for different widths L (measured 
in unit cells) and disorder strengths W (in units of the bulk 
gap A). Error bars are smaller than the symbol size, (a) and 
(b) Strong disorder: variation has flux-period hc/e and A is 
maximum at flux tp — n. Inset: Example of the exponential 
decay of the Green's function used to determine the local- 
ization length, (c) Weak disorder: Now, the hc/2e oscillation 
period is also apparent and localization lengths are very large 
(note log scale here). For wires much shorter than the zero 
flux localization length, the hc/e period will not be apparent, 
(d) Conductance Oscillation: Direct measure of conductance 
oscillation for a wire of length I — 100 and width L — 9, with 
strong disorder W = 2A. Maximum conductance occurs at 
7r-flux. 



even in the clean limit, metallic behavior appears at odd 
multiples of 7r flux, for a carefully tuned chemical poten- 
tial near the node. However, on raising the Fermi energy, 
the flux strength with the larger number of modes oscil- 
lates between flux zero and 7r. A robust response however 
is exposed by the presence of strong disorder, which we 
discuss next. 

Disordered System: We now consider adding disorder 
only on the surfaces of the model above, via a random 
onsite potential ViUi, (T invariant disorder), where rii 
is the charge density onsite, and Vi a random variable 
picked from a box distribution [-W, W]. We calculate 
the Greens function for a system composed of N layers 
in the z direction, and transverse size L x L, at energy 
E: GN(r,r';E). The position and spin coordinates are 
lumped into r. The transport properties of the surface 
states are characterized by extracting the localization 
length of the system in the quasi-one dimensional limit. 
Subsequently, we will directly calculate the conductiv- 



ity. The technical simplification for localization length A 
is that one only needs the Greens function between the 
first and last layers of the system, as the system length 
is increased. Denote by the matrix G^ N+ [, the Greens 
function between the first and last layers of a system with 
N + 1 layers, where the matrix indices refer to sites in 
the layer, and the energy label is suppressed. The local- 
ization length: A -1 = Urn ^log(\G\ M ,]| 2 ), which is 
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self averaging, is extracted by a linear fit (an example is 
shown in figure [2k inset) to the logarithm of the Greens 
function. The latter can be efficiently calculated [HI Q2] ■ 
Results are shown in figure [2] where parameters were 
chosen to obtain a bulk gap A = 2t, the hopping 
strength. System sizes with perimeter 4L, with L = 
6, 9, 12, 15 were studied and the quasi ID localization 
length X(<p) was extracted as a function of flux. We con- 
sider strong disorder, to obtain a localization length short 
enough to be measured: W — 2A for the first two and 
W = 4A for the last two sizes. The chemical potential 
was taken to be near the middle of the gap (E = 0), 
where the Dirac node appears in the clean limit. For 
these strong disorder strengths, results are nearly inde- 
pendent of chemical potential location inside the bulk 
gap. As seen in Figure [2^,,b a clear maximum in local- 
ization length is seen near ip = tt. Note, time reversal 
is already an approximate symmetry even for the small- 
est sizes, given the location of the maximum near ip = tt 
and the p — > 2tt — ip approximate symmetry. However at 
p = tt, localization length is finite, so time reversal sym- 
metry breaking enters here. Clearly, for larger widths 
the symmetry is more accurate, given the weaker fields 
prevailing on the surface states. If we denote: 

9<p = A(^)/4i 

this ranges from go = 1.6 to g^ — 6.3, as the flux is 
varied, for L — 15. Hence, for these parameters, wires 
with aspect ratios roughly within this range should ex- 
hibit conductance oscillations as a function of flux with 
a maximum at flux tt. This is explicitly checked below. 
We have also checked that the Zeeman splitting induced 
by the field for these sizes do not affect results qualita- 
tively, if their energy scale E z < 0.025A, which should 
be readily satisfied in experiments. 

Weak Disorder: In Figure [2H, we plot a system with 
weak disorder, where W = 0.5A, and the chemical poten- 
tial is tuned to E = 0.15 A, where there is a large density 
of states. The localization lengths are now significantly 
longer (note the log scale for A), and a prominent anti- 
localization feature is present near flux tp = 7r/2, Stt/2. 
The latter will contribute to an AAS hc/2e period in con- 
ductance oscillations. Since go — 125, at smaller aspect 
ratios the system is unaware that there is an even longer 
localization length at tt flux, and the most visible feature 
is likely to be the hc/2e AAS oscillations. This illus- 
trates the important role of strong disorder in observing 



the effect of interest. 

Conductance: Finally, in Figure [2ji, we present the 
conductance of a wire with length I — 100, cross sec- 
tion 9x9 and strong disorder strength W = 2 A. The 
conductance is extracted from the Greens function us- 
ing the Kubo formula, a = i^pTr [v z ImGv z ImG] with 
v z = -| [H, z] . The effect of leads is modeled by sand- 
wiching the disordered system between 4 x 10 3 layers of 
clean wire on either side. We take n — 0.15 A, to obtain 
a finite density of states in the leads. The localization 
length for these parameters is close to that in Figure [2J3, 
so I is intermediate between the localization lengths ob- 
tained there. A small imaginary part f ~ 5x 10~ 4 A is 
inserted in the energy to obtain finite results; results are 
insensitive to its precise value. Each data point is ob- 
tained by averaging log of conductance over 500 samples. 
Clearly, a conductance maximum at flux tt is observed. 
Note however, the hc/2e oscillations are more prominent 
here than in the localization length plots, and the over- 
all contrast in conductance is of order 0.1e 2 /h. This will 
increase for longer wires and wider cross sections, where 
the effects of T breaking at tt flux are less important, 
and localization of states away from this flux will set in. 
While a more extensive conductance analysis is left to 
the future, these results corroborate the basic picture. 

Scaling Analysis: How do these results scale to experi- 
mentally relevant system sizes, and in particular will the 
effects described here survive? As a reference, we note 
that a Bi 2 Se3 nanowires studied in Ref. [7J, was 350nm 
in circumference , and about 2/xm in length. While a 
direct conversion into unit cell lengths is difficult due to 
the extreme anisotropy of the material, we estimate as- 
pect ratios (length to circumference) between 5-10. The 
circumferences themselves are about 5-10 times larger 
in unit cell lengths, than those considered here. While 
in the truly two dimensional limit a metallic (symplec- 
tic metal) phase is expected, hence g[L — > oo) — > oo, 
this growth is slow. The scaling to larger widths is cap- 
tured by the beta function j3(g) — dlogg/dlogL, similar 
to the well known beta function for conductance. This 
function is known for spin-orbit metals in 2D [TJ]. For 
large g, the topological insulator surface displays identi- 
cal behavior [151 HI)] : i.e. 0(g) ~ 0.64/g. We estimate this 
to be reasonably accurate for g > 2 for TI surfaces, then 
g(rL) m g(L) + 0.64 log r. The localization length ratio 
for t = 10 times wider system compared to the L = 15 
case, will be g(L = 150) « 3.2 at zero flux , an aspect ra- 
tio that is easily exceeded. Note, the localization length 
at tt flux diverges more rapidly, since the effects of T 
breaking are weaker at larger widths. It is readily seen 
that this scales asymptotically as \(ip = tt) ~ L 4 (since 
the time reversal symmetry breaking strength scales as 
5 ~ 1/i 2 , and the localization length for weak T break- 
ing disorder scales as A ~ l/<5 2 )- Thus the oscillatory 
effects should remain visible at these longer scales, with 
strong disorder. 



Experimental Realization: In [7j, nanowires of Bi2Se3 
were subjected to a field along their length and an oscil- 
latory dependence of conductance 0(e 2 /h) with applied 
flux was reported. The oscillation period was consistent 
with surface states sensing a hc/e flux. However, the 
conductivity maxima were generally located at integer 
multiples of the flux quantum, rather than half integer 
multiples as predicted here. Therefore we believe a dif- 
ferent effect is at play there. The main ingredients re- 
quired to access the it flux effect discussed here is quasi 
one-dimensionality and strong disorder. The former can 
be achieved by studying narrower wires; eg. wires with 
half the width of those studied were produced in [7] . For 
the latter, one is aiming to reduce the metallicity of the 
surface layer, which in a Drude model is just (e 2 /h)k F l. 
This can be achieved either by decreasing the scattering 
length I, via greater surface disorder, or tuning the chem- 
ical potential to an energy with smaller carrier density. 

Conclusions: We described a topological property of 
the surface states of TIs. An Aharonov-Bohm phase of 7r 
strengthens the metallic nature of surface states, leading 
to a clear-cut transport signature. Conductance oscilla- 
tions with flux are expected in nanowires of topological 
insulator, with maxima at odd-integer multiple of ^o- 
The key requirements for observing this effect are quasi 
one-dimensionality and strong disorder, which we believe 
are achievable given current experimental capabilities. 

We thank H. Mathur, D. Carpentier, J. Moore, and 
G. Paulin for insightful discussions. This work was sup- 
ported by DOE grant DE-AC02-05CH11231. In indepen- 
dent work [18 , similar results are obtained using an ideal 
Dirac dispersion with twisted boundary conditions as a 
simple model of TI surface states experiencing magnetic 
flux. 



Supplementary Material 

1. Curved Surface Dirac Theory: Here we remind 
the reader of general formalism of defining curved 
spacetime Dirac theories, and then specializing to 
our case of interest where only the space is curved. 
We follow reference [17] . 

Consider flat 2+1 dimensional spacetime. We work 



9W 



with Euclidean time, so x 



it, and h 



c = 1. Then, the massless Dirac equation reads: 
(X)a=oT a ^ a ) ^ — 0- The matrices satisfy the an- 
ticommutation relations: {^ a , 7&} = 2S a b, where 
the right hand side is the Kronecker delta function, 
due to the Euclidean spacetime. The 7 matrices 
can just be taken as the Pauli matrices, and we can 
write this as a three vector 7 = (70, 71, 72). 

Now consider a curved spacetime, defined by coor- 
dinates x 11 metric g^, where the indices take values 
€ {0, 1, 2}. It is always possible to define a triad of 
vectors e M , with three components each, such that 



( • € v . One can define also a dual frame e M 
such that e M -e„ =8$, the Kronecker delta function. 
(This is like the relation between bases in real space 
and reciprocal space, for a crystal). Now, the gener- 
alization of the Gamma matrices can be defined as: 



ryfj- — e A" . ry ) or m component form 7^ = *)]„. 

The Dirac equation in curved space takes the form: 



e a la- 



7 M (^ + r M )* = o 

where d^ = d/d x n , and summation over repeated 
indices is assumed from here on. The T M matrices 
are the spin-connections, the main objects of inter- 
est here. They are given most compactly by the 
formula 



r„ = 



v /i7 " 
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(5) 


ui u + r^ x 




(6) 


9 va (dng a x + dxg,j, a - 


- d a g^ x ) 
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where the Eqn. [5] relates the spin connection to 
the covariant derivative V^. This is defined via 
the usual Christoffel symbols, T" x (Eqn. [7), which 
in turn are defined from the metric via the stan- 
dard formula (Eqn. u\. The spin connection can 
be shown to satisfy the equation: 



[1^,71 = ^+1^ 



(8) 



which indicates that it is responsible for parallel 
transport of spinors. 

We now specialize to the case when only space is 
curved. Then, we expect the metric tensor can be 
brought into the form g 00 = 1; g oi = 0, where the 
Latin indices i,j, k . . . will take on values G {1, 2}. 
This can be arranged if e° is unit modulus and or- 
thogonal to the other two tangent vectors. Note 
also, the only non- vanishing Christoffel symbols are 
ones with all Latin indices T l , k . Similarly the tem- 
poral spin connection vanishes r = 0. To read 
off the Hamiltonian, we write the Dirac equation 
as d^/dr = -if*, where H = 7o[7;(d; + I\-)]. 
We therefore identify (3 = 70, and using /3 2 = 1, 
7* = —if3a l . Substituting this into Eqn. |5j we 
obtain: 



r< = -\(mP + a k d t a k + T\ k a 3 a k ) 



(9) 



Actually, an even simpler form for the spin connec- 
tion is: 



L; 



1 



m/3 



(10) 



To derive this we first show (i) the spin connection 
is in the tangent plane, i.e. it can be expressed in 



terms of a, and then (ii) it is given by the expres- 
sion above. 

To do this, its useful to use the relationsjH] First 
consider the zeroth component 7 = /3; 

since all the relevant Christoffel symbols vanish. 
Now, if we can prove that I\ are in the tangent 
plane, then the commutator can be written as 
2f3Ti = —di(3 which gives us the desired result Eqn. 



10 We now prove the assertion that I\ lives in the 
tangent plane. To do this we write equations for 
atf. 

[T h a j } = d ia j +T\ k a k 

which can be derived from EqrJ8l and the definition 
a J = i(3~{ J . Now, if I\ only has components in the 
tangent plane, the right hand side is proportional 
to /3, i.e. it corresponds to the surface normal. We 
can check therefore that the right hand side has no 
components in the plane. Note, it corresponds to 
the vector c^e 5 + T J ik e k . However, we have men- 
tioned before that [c^e^],, = —T J ik e k . Therefore 
the only remaining component is perpendicular to 
the plane, as we required. 

To see this algebraically, we consider the inner 
product with a general in plane vector e;, which 
gives: 



e/ 



fte* 



1 a 



(11) 



5 k . This can be shown to 



where we used e fe • e; 
be zero. The critical input in this derivation is that 
diBj = djei, which is only true because these are 
tangent vectors, derived by differentiating the three 
dimensional coordinates of the surface e* = dr/dx l . 

We note identical results are obtained if we choose 
the basis vectors as a fixed linear combination of 
tangent vectors: i.e. if e^ = cos 6e\ + sin 9e2 and 
e 2 = — sinflei +cos#e2. In particular, a natural 
choice is 9 — tt/2, which gives a Dirac Hamiltonian 
in flat space of the form H = —ivph(a y d x — a x d y ). 

We now discuss a general argument which fixes the 
half integer offset of momentum quantization, and 
goes beyond the particular Dirac model chosen. 
Note, by time reversal invariance fixes this offset 
to be 7r or zero. In the latter case, a gapless one 
dimensional mode is present that is protected by 
time reversal symmetry. Now consider a slightly 
different geometry, an annular cylinder with radii 
i?i, i?2- Take R2 — > 00 so we now have a hollow 
inner cylinder. If we now shrink the inner radius 
R% —>■ 0, one just obtains the bulk topological in- 
sulator. Hence, we must have that the modes on 



the inner cylinder were gapped, which implies 7r 
phase shift. Note, although the cylinder is 'flat' 
in the sense it has vanishing Gaussian curvature, a 
nonzero spin connection leads to the Berry's phase 
of 7r. It is in this sense that we refer to this phase as 
being surface curvature induced - strictly it refers to 
a topological property of the surface (i.e. the possi- 
bility of looping around the cylinder while avoiding 
the inserted magnetic field). 



2. Model: We use the Fu-Kane-Mele model of topo- 
logical insulators on a diamond lattice Eq. [4] where 
the d ik are nearest neighbor unit vectors connect- 
ing a pair of second neighbor sites ik. We choose 
nearest neighbor hooping to be strong along direc- 
tion (111) with strength t\ = It and the remaining 
three bonds to be of equal strength t. The spin 
orbit interaction is taken to be Xso = 0.25£. The 
'cylinders' used for the computations are actually 
parallelepipeds, with the long axis z being along 
(110), and the cross section axes being along (Oil) 
and (101). Note, the single Dirac node at the sur- 
face of the long faces is located at k z — n. This 
does not affect any of the results since this wavevec- 
tor is along the propagation direction, and can be 
ignored. Note, symmetry of the diamond lattice 
relates the surface Dirac nodes along the two dis- 
tinct surfaces, hence they are at the same energy. A 
more general situation is when there is no particular 
relation between the two - however, the topologi- 
cal properties should however remain unchanged, 
including the anomalous Aharonov-Bohm oscilla- 
tions. 



3. Green's Function Calculation: The matrix 
G\ N+ [(E), the Greens function between the first 
and last layers of a system with N + 1 layers at 
energy E, where the matrix indices refer to sites 
in the layer. This can be efficiently calculated fT3"] . 
assuming only nearest neighbor hopping between 
layers, denoted by the matrix t and t^. If the 
single layer Hamiltonian is Hjv+i, then the Greens 
function within this layer when attached to the 
remaining N layers is just: 



G 



(N+l) 
N+1,N+1 

n (AT+l) 



(£7 + w/)-HAr +1 -ttG^t 



Vj l,JV b,J Af+l,JV+l 



(12) 
(13) 



where the energy label is suppressed. The last 
equation gives us the desired Greens function. 
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